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We study the cosmic microwave background (CMB) anisotropy due to spherically symmetric 
nonlinear structures in flat universes with dust and a cosmological constant. By modeling a time- 
evolving spherical compensated void/lump by Lemaitre-Tolman-Bondi spacetimes, we numerically 
solve the null geodesic equations with the Einstein equations. We find that a nonlinear void redshifts 
the CMB photons that pass through it regardless of the distance to it. In contrast, a nonlinear lump 
blueshifts (or redshifts) the CMB photons if it is located near (or sufficiently far from) us. The 
present analysis comprehensively covers previous works based on a thin-shell approximation and a 
|. linear/second order perturbation method and the effects of shell thickness and full nonlinearity. Our 

D . results indicate that, if quasi-linear and large (> lOOMpc) voids/lumps would exist, they could be 

observed as cold or hot spots with temperature variance > 10 _5 K in the CMB sky. 
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■ PACS numbers: 98.80.-k, 98.70.Vc, 04.25.Nx 

^ ; I. INTRODUCTION 

Recently, much attention has been paid to generation of the cosmic microwave background (CMB) anisotropy 
due to nonlinear evolution of the gravitational potential, which is called the Rees-Sciama (RS) effect [l|. It has 
| been argued that the RS effect only affects the angular power spectrum of the CMB anisotropy at relatively 
£NJ . small angular scales I > 3000 0, 0, El- However, recent discoveries of the CMB anomalies such as octopole 
^ ' planarity, the alignment between quadrupolc and octopole components [f| , anomalously cold spots on angular 
', scales - 10° @], and asymmetry in the large-angle power between opposite hemispheres Q hint that the RS 
effect due to large-scale structures could affect the CMB anisotropy at large angular scales as well Q. This 
\ possibility is also indicated by a recent report that the density of extragalactic radio sources as projected on 
• ■ the sky is anomalously low in the direction towards the cold spot in the CMB map @. 

The signatures of the RS effect due to nonlinear voids/lumps in the Friedmann- Robertson- Walker (FRW) 
universe without a cosmological constant has been extensively studied in the literature 0, d, 01 • Recently, the 
RS effect due to a quasi-linear void/lump in the FRW universe with a cosmological constant has been studied 
using a thin-shell approximation (Toj and a second-order perturbation method pd| . In order to check the validity 
. £h ' an< ^ cons i s t enc y; h is of great importance to extend the analyses to solve the Einstein equations without relying 
k> , on these approximations. 

■ In this paper, we study the RS effect due to nonlinear structures in flat universes with a cosmological con- 
Ci ' stant by solving the Einstein equations, which incorporate the fully nonlinear regime. Specifically we model a 

compensated spherical nonlinear void/lump by a family of Lemaitre-Tolman-Bondi (LTB) spacetimes, and nu- 
merically solve the null geodesic equations with the Einstein equations. In §2, we derive the Einstein equations 
and null geodesic equations for LTB spacetimes and model a compensating spherical void/lump with a smooth 
mass density profile. In §3, we show some numerical results. §4 is dedicated to concluding remarks. 
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II. MODEL AND BASIC EQUATIONS 
A. Lemaitre-Tolman-Bondi spacetime 

We consider a spherically symmetric spacetime with dust and a cosmological constant A, which satisfies 
Einstein equations, 

G M „ + Kg^ u = BirGpu^Uy (2.1) 

where g llUl G^,,, G, p, and arc the Ricmannian metric tensor, the Einstein tensor, the gravitational constant, 
matter density, and the fluid 4-velocity, respectively. 

In spherical coordinates (t, r, 9, </>), the general solutions are represented by Lemaitre-Tolman-Bondi (LTB) 
metric, 

ds 2 = -dt 2 + — ^\ dr 2 + R 2 (t, r)(d9 2 + sin 2 6d<p 2 ), (2.2) 
1 + f{r) 



which satisfies 



2Gm{r) A 



R=— ^ + -R 2 + f(r), (2.3) 
m'(r) 

where ' = d/dr and ' = d/dt. The solutions contain two arbitrary functions, /(r) and m(r). If we we give 
'A, p(ti,r), and the local Hubble parameter H(ti,r) = R(ti , r) / R(ti , r) at the initial time t = ti, m(r) and 
f(r) are determined by (|2.3[) and (|2.4[) . The radial coordinate r has a gauge degree of freedom, r — * r' = [any 
function of r]; here we define r as the areal radius at the initial time: R(ti, r) = r. 

Once m(r) and f(r) are determined, the evolution of R is given by (|2.3|) numerically. Differentiating (|2.3p 
with respect to r, we obtain 

* = _L(^_-* + f + ^), 

which is the evolution equation of R' . Although R' can be calculated also by the finite difference of R with 
respect with r, the integration of f !2.5[) with respect to t gives better precision for R' . 



B. Modeling a Void/Lump 



Our model is composed of three regions: the outer flat FRW spacetime (V+), the inner negatively/positively 
curved FRW spacetime (V_), and the intermediate shell region (V s ). In V±, the field equations (|2.3[) and (|2.4p 
reduce to the Fricdmann equations, 

(2.6) 
(2.7) 
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where G = yj is a constant. Here r = r± denotes the boundary between V± and 14 whereas o±, ff±, 
and p± simply mean the quantities in V± . 

The shell V s is constructed by the LTB spacetime in such a way that m(r) and f(r) are continuous through 
V±. At the initial time t = ti, we assume that p_ w p + , _ff (tj, ?*)=const., and the matter density profile is given 
by 



P- 

Pc 



p(r) = < 



— — (3Xi - lOATi + 15X_ 
16 

Pc -{3Xl - 10X1 + 15X+ 
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for r < r_ , 

p_ for r_ < r < r c , 

p c for r c < r < r + , 

for r > r+, 



(2.8) 
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FIG. 1: Examples of initial and evolved profiles of p(t,r). (a) and (b) correspond a void (S < 0) and a lump (5 > 0), 
respectively. 



with ry = 



r + + r_ 



r + — r_ 



X+ = 



_ T — r c =p w/2 



w/2 



(2.9) 



Among the parameters above, p c = p(r c ) cannot be fixed in advance. It is determined as an eigenvalue of 
Einstein equations (|2.3[) and (|2 .4[) with the boundary condition m(r + ) = 47rG / o+(a+r + ) 3 /3 and /(r+) = 0. We 
define the times t\, t%, and t$ as follows: the photon passes r = r + from V+ to V s at t = t\, passes r = r_ 
from V" s to V- at i = t^, passes the other side of r = r_ from V— to at i = t3 and passes the other side of 
r = r + from V s to V+ at i = t±. Examples of initial and evolved configurations of p(t, r) are shown in Fig. 1. 

Out model parameters are the background density parameter, the density contrast, the physical radius in 
unit of the Hubble radius, and the width of the shell in unit of the comoving radius of the void/lump, 



_ 8nGp + (t 4 ) _ P-{U) , ul , \ jj / j \ , 



(2.10) 



at the exit time t 4 of the photon. The initial parameters O^, Si, and r c are obtained by iterative integration of 
the field equation (|2~3|) with IpTTOj) , 



C. Temperature anisotropy 

The wave 4-vector fc M of a photon satisfies the null geodesic equations, 

dx 1 * 

h = , k^k^ = 0, 



dk^ 

^ + F^=0, 



(2.11) 
(2.12) 



where A is an affine parameter. In what follows, we only consider a CMB photon which passes the void/lump 
center, r = 0. Then the geodesic equations (|2.1ip and (|2.12|) with the metric (|2.2[) yield 



— - k l — - v k 9 - kf - n 
dx - ' dx ~ ' _ " u ' 



■^^k*, e = sign| 
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(2.13) 
(2.14) 

(2.15) 
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For the period t% < t < £3 in V— , the evolution of k l _ and the crossing time arc given by 



1 f ta dt 2 



k*_ oc — , / — = — arcsinh(Cr-). (2.16) 



For the periods t\ < t < t% and £3 < t < t& in V s wc numerically solve the geodesic equations and the field 
equations simultaneously. First, we discretize the radial coordinate into N elements, 

n = r- + (i - l)Ar, i = l,...,JV, Ar = ^— (2.17) 

and R(t,r) into i?i(i) = R(t,ri). 

Next, we rewrite the geodesic equations (|2.13[) - (|2 . 15[) and the field equations (|2.3[) and (|2.5[) as differential 
equations of r, 



eR' 



dr VT + 7 
dfc* ei?' 



(2.18) 

fc*, (2.19) 



*( 9rr kn = ^*r, (2.20) 

dr 2 

where fi, i?' and [ ^ ) are given by (12.3H . (12.51) and (12.181) . respectively. 
\ dr J 

Finally, we carry out numerical integration of (12.18[) . (|2.19p . (|2.2ip . and (I2.22p using the fourth-order Rungc- 
Kutta method to obtain the solutions of t(r), fc*(r), Ri(t(r)), and i?^(f(r)). To estimate the numerical precision, 
we also numerically solve Eq. (|2.20j) and check how the solution satisfies (violates) the null condition (|2.14|) . 

The energy of a photon passing through the homogeneous background without a void/lump is 

fc* oc — . (2.23) 
a+ 

Then the temperature fluctuation caused by a void/lump can be written as 

— = £ - 1. (2-24) 

T k\ V ' 



III. RESULTS 

Figure 2 shows temperature fluctuations of photons passing through a void/lump for comoving observers at 
each r = constant. The amplitude of fluctuations temporarily increases to \AT/T\ ~ 10~ 3 , but it finally reduces 
to ~ 10 -5 at the edge of the shell because the mass of the void/lump is compensated. 

In what follows, we discuss only the values of AT/T measured by a comoving observer outside a void/lump. 
For a void, as Fig. 3(a) indicates, the temperature fluctuation AT/T is always negative regardless of the values 
of H4. For a fixed radius, \ AT/T\ decreases as the width w/r c of the void shell increases. Wc find that our results 
are consistent with those for a thin-shell homogeneous void in the quasi- linear regime (Io| . To see nonlinear 
effects, in Fig. 3(b)(c) we plot AT/T obtained from a linear perturbation analysis, a second order perturbation 
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FIG. 2: Temperature fluctuations of photons passing through the center of a large void (a) and a large lump (b) for 
comoving observers at each r — constant. The parameters are £4 = TO-3, Q.4, = 0.24, Ri{r c ) = Q.IH^ 1 , and w/r c = 0.1. 
The subscript 4 denotes quantities at the time £4 when a CMB photon exits the edge of a void/lump. The arrow indicates 
the traveling direction of a CMB photon. 




FIG. 3: Temperature fluctuations of CMB photons passing through the center of a large void with Ri(r c ) = 0.1ff 4 1 for 
a comoving observer outside the void, (a) shows AT/T as a function of fk for 84, — —0.3. The dotted line indicated by 
"thin shell" shows AT/T for the thin-shell model [Tfl]- (b) and (c) show AT/T as a function of -64 for Q 4 = 0.24 and 
for Q4 = 0.9, respectively; we put w/r c — 0.3 for both cases. The dotted lines and the dashed lines represent the values 
obtained from a linear perturbation analysis, and a second order perturbation analysis [TJ], respectively. 
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FIG. 4: Temperature fluctuations of CMB photons passing through the center of a large lump with Ri(r c ) = O.IH^ 1 
and w/r c = 0.15 for a comoving observer outside the lump. The dotted lines and the dashed lines represent the values 
obtained from a linear perturbation analysis, and a second order perturbation analysis [ll|], respectively, (a) shows AT/T 
versus Q4 for S4 = 0.3. (b) and (c) show AT/T versus £4 for Q.4 — 0.24 and for Q4 = 0.9, respectively. 



analysis, and our fully nonlinear analysis. We find that higher-order effects are important and they enhance 
I AT/T I for a void, particularly for large 

For a lump, as Fig. 4(a) indicates, the temperature fluctuation AT/T is positive for low background density, 
but it can be negative for high background density. In other words, lumps at low-z blueshift the CMB photons, 
whereas lumps at high-z redshift them. This behavior of AT/T in the quasi- linear regime can be interpreted 
as follows. First, let us consider a perturbative case \8i\ <C 1 for which the linear integrated Sachs- Wolfe (ISW) 
effect [lH dominates AT/T. Then one can show that AT/T vanishes for ^4 = (de Sitter) because no matter 
fluctuation exists, and for = 1 (Einstein-de Sitter) because the Newtonian gravitational potential freezes in 
the Einstein-de Sitter universe. Therefore, AT/T cannot be a monotonic function of ^4 for < O4 < 1. In 
fact, as one can see in Fig. 4(a), the ISW contribution has a peak as a function of ^4 for a fixed 64. Next, let 
us consider a quasi-linear case 0.1 < < 1. For small Q4, the nonlinear RS effect is not important because 
matter fluctuations are small. However, for large ^4, the nonlinear RS effect dominates the linear ISW effect, 
which vanishes for = 1. Our numerical analysis shows that the nonlinear RS effect reduces the temperature 
of the CMB photons, which reconfirms the previous semi-analytic result for lumps in the Einstein-de Sitter 
universe Thus, one can interpret that negative AT/T in the Q4 = 1 background for a void/lump (in Fig. 
3(a)/4(a)) is caused by the nonlinear RS effect alone. It should also be noted that our result is consistent with 
the previous one obtained from a second order perturbation analysis for a void/lump in accelerating universes 

a 

Fig. 4(b) (c) shows nonlinear effects for a lump: higher-order effects are still important, but they reduce the 
amplitude of AT/T in contrast to the case for a void in Fig. 3(b)(c). 
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IV. CONCLUDING REMARKS 

We have studied the CMB anisotropy caused by spherically symmetric nonlinear structures in flat universes 
with dust and cosmological constant. Specifically, by modeling a time-evolving spherical compensated void/lump 
by Lemaitre-Tolman-Bondi spacetimes. we have solved the null geodesic equations with the Einstein equations 
numerically. 

We have found that a nonlinear void redshifts the CMB photons that pass through it regardless of its location. 
In contrast, a compensated nonlinear lump blucshifts (or redshifts) the CMB photons if it is located near (or 
sufficiently far from) us. 

Our result for the temperature anisotropy due to a void is roughly consistent with the previous one based on 
a thin-shell approximation. We have also shown that \AT/T\ decreases as the shell thickness increases for fixed 
6. 

We have checked that our results are also consistent with the ones based on a linear/second order perturbation 
method for small \d\. It turned out that nonlinear (higher-order) effects are important even in the quasi-linear 
regime |<5| > 0.3. 

Our results indicate that, if a quasi-linear (\5\ ~ 0.3) and large size (R ~ O.IH^ 1 ) void/lump could exist, 
they would be observed as a cold or hot spot at the level of AT/T ~ 10~ 5 in the CMB sky. In such a case fully 
nonlinear and relativistic analysis is necessary. 
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